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The lightest supersymmetric particle, most likely the neutralino, might account for a large fraction 
of dark matter in the Universe. We show that the primordial spectrum of density fluctuations in 
neutralino cold dark matter (CDM) has a sharp cut-off due to two damping mechanisms: collisional 
damping during the kinetic decoupling of the neutralinos at about 30 MeV (for typical neutralino 
and sfermion masses) and free streaming after last scattering of neutralinos. The last scattering 
temperature is lower than the kinetic decoupling temperature by one order of magnitude. The cut- 
off in the primordial spectrum defines a minimal mass for CDM objects in hierarchical structure 
formation. For typical neutralino and sfermion masses the first gravitationally bound neutralino 
clouds have to have masses above 1O _7 M0. 

PACS numbers: 14.80.Ly, 98.35. Ce, 98.80.-k, 98.80.Cq 



I. INTRODUCTION 

Recent observations of the cosmic microwave back- 
ground (CMB) anisotropies are consistent with a key 
prediction of inflationary cosmology: The universe ap- 
pears to be spatially flat Q. However, only a small per- 
centage of the mass that is needed to account for the 
critical energy density of the universe comes in the form 
of baryons. Recent observations of primordial deuterium 
and other light elements suggest that the baryonic mass 
density is w b = 0.019 ± 0.0024 §, which implies that 
only about 5 per cent of the mass in the universe is bary- 
onic. The remaining mass is assumed to be a mixture 
of different forms of yet unknown dark matter and dark 
energy. However, we do have evidence, mainly from the 
study of large scale structures, about the properties of 
dark matter. 

Cold dark matter (CDM) by definition has a non- 
relativistic equation of state at the beginning of struc- 
ture formation around the matter-radiation equality || . 
For successful structure formation an important fraction 
of the dark mass has to be cold dark matter. Although 
purely baryonic matter and hot dark matter (relativistic 
equation of state at matter-radiation equality) models 
have been ruled out long ago Q], a model with a cos- 
mological constant and baryonic matter only provides a 
good fit to the recent CMB observations (|] . When com- 
bined with other cosmological observations it turns out 
that the small sound speed (at photon decoupling) of 
the baryonic matter can only be compatible with the ob- 
served multipole moments if the universe is closed ||. 
Moreover, this model does not provide enough power at 
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small scales to explain the observed distribution of galax- 
ics §. 

The most important feature of CDM is hierarchical 
structure formation, i. e., small structures form first and 
grow to larger structures later. A natural question ob- 
viously is: What is the smallest mass scale [@ of CDM 
structures? 

Since the nature of CDM is unknown, the answer to 
this question will not be the same for different CDM 
candidates. Natural candidates are particles that are 
predicted by extensions of the standard model of parti- 
cle physics. A CDM candidate should be (meta-)stable, 
(electrical and color) neutral, and heavy.The minimal su- 
persymmetric standard model (MSSM) with the assump- 
tion of conserved R-parity provides an excellent candi- 
date: the lightest neutralino x®, which probably is the 
lightest supersymmetric particle (LSP) 

The time of kinetic decoupling of CDM depends on the 
nature of CDM During kinetic decoupling colli- 

sional damping is the dominant mechanism. Once CDM 
is fully decoupled from the radiation fluid, damping due 
to free streaming happens. Interesting general consider- 
ations on damping mechanisms for CDM have been pub- 
lished by Boehm, Fayet, and Schaeffer [|llj recently. 

After the neutralinos decouple chemically (at about 
T c( j ~ M-j/20) they remain in kinetic equilibrium due 
to frequent scattering with particles from the radia- 
tion fluid. After the QCD transition (at ~ 160 MeV) 
neutralino-lepton scattering is the most important pro- 
cess. The neutralinos decouple kinetically once the relax- 
ation time t becomes comparable with the Hubble time 
in = H , which happens, depending on the parameters 
of the MSSM between 10 MeV and a few 100 MeV. Once 
collisions of neutralinos with particles from the radiation 
fluid cease, the equation of state becomes nonrelativistic 
(P 0) and neutralino matter starts its life as cold dark 
matter ||^||. 



1 



In the present work we calculate the temperatures of 
kinetic decoupling and last scattering of neutralino CDM 
for the case of a bino-like neutralino. A first estimate 
of the kinetic decoupling temperature, based on dimen- 
sional arguments, was given by Schmid, Schwarz, and 
Widerin ({J , which has been confirmed recently by more 
detailed calculations by Chen, Kamionkowski, and Zhang 
|jl2"|| . However, the authors of Ref. (IJ] ignored the fact 
that the relevant time scale that has to be compared 
to the Hubble time is the relaxation time, rather than 
the collision time. It was shown explicitly that photon- 
neutralino scattering is suppressed by several orders of 
magnitudes, compared to lepton-neutralino scattering. 

During the process of kinetic decoupling, collisional 
damping can smear out primordial fluctuations in neu- 
tralino CDM below some mass scale M<±- Free streaming 
gives rise to additional damping below Mf s (t), which de- 
pends, in contrast to Md, on time. Both damping mech- 
anisms together give rise to a sharp cut-off in the primor- 
dial power spectrum of neutralino CDM, that typically 
lies at M ~ 10~ 7 M Q at the time of matter-radiation 
equality. We have presented preliminary estimates in 
Ref. |jig. In Ref. @ it was pointed out that the estimate 
of induced damping found in [ fi"l| , is wrong by several or- 
ders of magnitude mainly because the cross section for 
elastic scatterings of photons with neutralinos has been 
overestimated. 

We also show that bulk viscosity, besides shear viscos- 
ity, can not be neglected (as has been done in Jllfl ) in 
the situation when a nonrelativistic component decou- 
ples from a radiation fluid. At first sight this is a sur- 
prising result since bulk viscosity usually goes along with 
the transfer of energy to internal degrees of freedom or 
with particle production. None of these mechanisms is 
available here. However, CDM and radiation have to be 
treated as two separate fluids and the bulk viscosity of 
the CDM fluid just reflects the energy dissipation from 
the CDM fluid to the radiation fluid, which is however 
a negligible effect for the radiation fluid since the energy 
density of the CDM fluid is tiny compared to the energy 
density of the radiation fluid at kinetic decoupling. On 
the other hand the heat conduction (which has been con- 
sidered in jll]]) can be neglected for the CDM fluid. The 
reason is simple, the neutralinos are too slow. 

The paper is organized as follows: A short summary of 
mass limits and our assumptions about the lightest neu- 
tralino is given in Sec. II. Then we review the simplest 
calculation of chemical decoupling for pedagogical rea- 
sons and compare that with our detailed calculation of 
the kinetic decoupling and last scattering temperatures 
(Sec. III). In the following section (IV) we introduce 
CDM as an imperfect fluid, along the lines described in 
Jl3| p"5| . The kinetic theory for the description of CDM 
is explained in Sec. V, and the coefficients of transport 
are calculated in Sec. VI. For this purpose we general- 
ize the program by Weinberg [ fl5| and Straumann |l6[ to 
the situation of a nonrelativistic component that decou- 
ples from a relativistic fluid [Weinberg and Straumann 



treat the problem of decoupling of a relativistic com- 
ponent (photons) from a nonrelativistic fluid (baryons)]. 
This finally allows us to calculate the damping scale from 
kinetic decoupling (Sec. VII) and free streaming (Sec. 
VIII). We conclude with a short discussion of the impli- 
cations of our findings. The relevant cross sections are 
calculated in App. A and some useful thermodynamical 
relations can be found in App. B. 

II. NEUTRALINOS 

A direct lower limit for the neutralino mass is pro- 
vided by the LEP experiments, > 37 GeV for any 
tan (3 and sfermion mass ||l7|| . Reasonable assumptions 
(universal soft supersymmetry-breaking scalar masses 
at some higher scale) within the MSSM (constrained 
MSSM) and taking results of the Higgs searches into 
account raises the lower limit to about 50 GeV |]l7j , |l8|| . 
Incorporating also constraints from b — > sj decays and 
assuming that neutralino dark matter is cosmologically 
interesting (0.1 < ui x < 0.3) a lower limit as high as 
Mx — 140GeV can be derived [Q. The cosmological 
upper limit also gives rise to an upper limit on the neu- 
tralino mass. It is essential to include the effects of the 
next-to- lightest supersymmetric particles and coannihila- 
tions, as well as the contribution from poles and thresh- 
olds properly |l8|-po|| . A detailed analysis gives M x < 600 
GeV jl8). Since many untested, although reasonable, as- 
sumptions go into these limits, we decided for this work 
to assume M x > 50 GeV. 

The neutralinos are linear combinations of the neutral 
gauginos and the two higgsinos of the theory, i. e., 

Xl = Z n B° + Z 12 W° + Z 13 H° + Z U H° (1) 

expressed in terms of mass eigenstates. The Zij , j S 
{1,2,3,4} are elements of a real orthogonal matrix which 
diagonalizes the neutralino mass matrix. In most of the 
parameter space of the constrained MSSM the LSP is a 
B°. We assume Xi = X ~ B°. 

For a pure bino the interaction with a standard model 
fermion F is given via the exchange of the related left- 
or right-handed sfermion Flu a s follows: 

C FF - = -V2g F {b F F L V R - c F F r Vl) + h.c. , (2) 

where g is the electroweak coupling constant, Vl,r de- 
notes the left and right chiral projection operator. The 
left and right chiral vertices are given by 

Y F 

b F = Zn— tan6> w + Zi 2 T 3F , (3) 
cp = ZuQptanOyj . (4) 

Here Yf, T 3 f and Qf are the weak hypercharge, isospin 
and electrical charge of the involved fermions. 
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III. CHEMICAL AND KINETIC DECOUPLING 

There is a large difference between the temperature of 
chemical decoupling, T cc j, and the temperature of kinetic 
decoupling, X^d, of neutralino cold dark matter. This is a 
characteristic feature of weakly interacting massive parti- 
cles (WIMPs) . Chemical decoupling (freeze-out) fixes the 
relic abundance of neutralinos and therefore the present 
value of ricDM^i 2 - Before kinetic decoupling the neu- 
tralinos are tightly coupled to radiation, after kinetic de- 
coupling the neutralinos acquire the properties of CDM. 
Namely, the neutralinos interact with radiation only via 
gravity and their pressure is negligible compared to their 
energy density well before matter-radiation equality. 

Let us first review the process of chemical decoupling, 
which is a useful warm up for the kinetic decoupling that 
is explained subsequently. We assume that the neutralino 
is the bino, which reduces the number of free parame- 
ters to the bino mass, M% and to the universal sfermion 
mass, Mp. For a more complete picture including hig- 
gsino admixture, thresholds, poles, and coannihilations 
see HHg. 

At T > T c d neutralinos are kept in chemical equi- 
librium with all standard model fermions F in the heat 
bath at temperature T via annihilation processes x+x 
F + F. From (^) one can calculate the annihilation rate 
iov X + x^F + F @, 



annihilation rate r ann becomes comparable to the ex- 
pansion rate H of the universe neutralinos no longer 
find other neutralinos to annihilate. We use the con- 
dition r ann = H to define the temperature of chemical 
decoupling T c d. Solving this equation iteratively yields 
(x ee Mx/T) 
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as long as the bino mass is well below the top mass, but 
large compared to the bottom mass. In deriving (f?j) we 
assumed equal masses for all sfermions. Exploring the 
parameter space of the MSSM we typically find x c( i « 25, 
cf. figure [l]. 
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1. The chemical decoupling as a function of the 



Here, (. . .) denotes thermal averaging and v is the 
Moeller-velocity. In order to obtain we expanded 
ii(j ann for small mp/M^ and small v. More details can 
be found in appendix A. Note that the first term in 
the square brackets contributes practically only for top 
quarks (F — t). However, we will assume below that 
< mt such that the second term will be dominant 
in our estimate. We neglect annihilation of neutralinos 
into final states containing gauge and Higgs bosons, such 
as xx <-> {WW, ZZ, HH, HW, HZ}, since these channels 
are particularly important for Higgsino-like and mixed- 
state neutralinos, but are subdominant when compared 
to the fermion-antifermion channels in the case that the 
neutralino is mostly a gaugino [ fL9[ . Since we restrict our 
attention to a pure bino, there is no contribution from 
diagrams with Z a exchange at tree level and therefore 
the Z° pole does not invalidate our estimate below. 

As the universe expands the temperature eventually 
falls below the neutralino mass and the number den- 
sity riy. of neutralinos decreases exponentially. Once the r c i = \_. ( v &c\(El)) (T) ^l(T) 



sfermion mass Mp for three values of the neutralino mass 
Mjj = 50, 100, 150 GeV (increasing from bottom to top). 

The relic abundance of neutralinos is now easily ob- 
tained as n-g(To) = n^(T c d)s(To)/s(T c( i), where To = 
2.725 K and s denotes the entropy density of the universe. 
It is a good approximation to use the equilibrium distri- 
bution for the number density at T c d, although in a more 
advanced treatment the corresponding kinetic equation 
should be solved. From the number density n^(T ) we 
may easily compute oj^ = fl^h 2 , which is plotted in Fig. ^ 
as a function of for typical values of the sfermion 
mass. 

Below T c d the neutralinos are kept in local thermal 
equilibrium via elastic scattering processes x + F ~* X + 
F. After the QCD phase transition only leptons L remain 
as scattering partners for the neutralinos. We neglect 
scatterings with pions, which is important for T > m t 
only. It will turn out that in most cases T^d <C m^. 
Scattering with nucleons is not important due to the tiny 
number density of baryons. From (^|) one can calculate 
the rate of elastic scatterings x + L —> x + L j2^]. We 
find 
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M L 



T 2 n L (T). (9) 



El denotes the energy and til the number density of the 
lcptons. In deriving (0) we approximate the Mandelstam 
variable s w M^ 2 + 2M X E L . Note that the Moeller ve- 
locity in this case is v ~ 1 to a very good approximation. 




FIG. 2. The relic abundance of neutralinos expressed by 
lu x = ft 2 as a function of the neutralino mass for differ- 
ent values of the sfermion mass Mp = 150, 200, 250, 300, 400 
GeV. The sfermion mass increases from the bottom to the 
top. The dark shaded region is excluded by the conservative 
assumptions: Q < 1 and h < 0.8. The light shaded region 
indicates typical values of uj in a ACDM model. 

In analogy to the chemical decoupling the condition 
r o i = H defines the temperature at the time the last 
elastic interactions between neutralinos and the rest of 
the Universe take place. This last scattering temperature 
is given by 



The leptons are kept in local thermal equilibrium through 
the frequent interactions among themselves and the 
equipartition theorem gives El = 3/2 T. Comparing 
the rms momentum transfer with the typical neutralino 
momentum p x we find Apx/Px = \/3/2 T/M^ <C 1. 
This means that a huge number N(T) of elastic scatter- 
ings is needed to keep or to establish thermal equilibrium, 
N(T) = p x /Ap x = \fW* M x /T. We can now estimate 
the relaxation time as 



r(T) 



Mx 
T 



Tcoll 



(12) 



Note that r(T) - 1/T 6 . 

The kinetic decoupling of the neutralinos happens 
when the relaxation time r becomes comparable to the 
Hubble time 1/H. We denote the corresponding temper- 
ature by Tkd, which is given by 



Tkd = 



1.2 x 10" 
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M X (M? - M|) 2 



-1/4 
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where we assumed that all leptons except the tau are 
relativistic, but we neglected the contribution of pions 
which arc important at temperatures of about 130 MeV. 
Above the QCD phase transition at about 160 MeV much 
more interaction partners are available and our formula 
should be modified. Exploring the parameter space of the 
MSSM we typically find that T kd is of the order 10 MeV 
to 100 MeV, cf. figure |. For M x = 100(150) GeV and 
M £ = 200(250) GeV we find T kd = 28(36) MeV, whereas 
the chemical decoupling for the same set of parameters 
happens at T c( j = 4.0(5.9) GeV, a difference of more than 
two orders of magnitude. 
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Typical values are 1 MeV to 10 MeV, e.g., T ls = 2.3(2.5) 
MeV for M x = 100(150) GeV and M~ L = 200(250) GeV. 
However, this is not the temperature at which neutrali- 
nos decouple kinetically. The kinetic decoupling temper- 
ature is defined through the relaxation time r, rather 
than by the collision time r co ii = l/r c i. This can be 
easily understood by the following argument. 

The relaxation time r, i.e., the time neutralinos need 
to return to local thermal equilibrium after a deviation 
from it, can be estimated from the typical number of 
scatterings that is needed to change the momentum of the 
neutralino significantly. The typical momentum transfer 
in a single elastic scattering event is tiny compared to 
the average momentum of the neutralinos. This is easily 
seen from the averaged Mandelstam variable t, 



1 



dcr c 
dt 



tdt 



2 El 



(11) 



> 



200 
175 
150 
125 



■a 100 

H 75 
50 
25 




200 



400 



800 



1000 



Ml 



600 

[GeV] 

FIG. 3. The temperature of kinetic decoupling of neutrali- 
nos from radiation as a function of the sfermion mass for 
M x = 50, 100, 200 GeV (bottom to top). 

The large difference between T c d and Tkd is mainly 
due to the different target densities in the annihilation 
[Eq. (g)] and elastic scattering rates [Eq. (g)]. For annihi- 
lations the target density is given by the number density 
ri^ of neutralinos in the universe. The number density 
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of neutralinos is suppressed by the Boltzmann factor at 
chemical decoupling. In contrast the target density for 
elastic scattering processes is given by the number den- 
sity of all relativistic leptons. 

During the cooling from P c d to P^d the leptons behave 
as a perfect radiation fluid which tries to keep the neu- 
tralinos in thermal equilibrium through clastic scattering 
processes. The neutralinos on the other hand may be de- 
scribed as a nonrelativistic, imperfect fluid. 



IV. CDM AS A FLUID 

For temperatures T > T c( j the particle content of the 
universe may be described by a single radiation fluid 
which is in local thermal equilibrium. For temperatures 
Tkd < P < T c d the radiation fluid is tightly coupled to 
the CDM fluid. Hence, both fluids have the same tem- 
perature and velocity four-vectors. Around p^d the CDM 
fluid starts to decouple kinetically from the radiation 
fluid and becomes an imperfect fluid. The departure from 
local thermal equilibrium is generated by dissipation, i.e. 
by shear and bulk viscosity (we show below that the co- 
efficient of heat conduction vanishes). For temperatures 
T < Tkd both fluids are decoupled and the CDM fluid 
is freely streaming. Since Ocdm = (a/a cq ) fl ra d <C fi ra d 
for P 3> P eq , the radiation fluid remains in local thermal 
equilibrium throughout the decoupling process. 

The current density and the energy-momentum tensor 
of the radiation fluid (R) are given by 



J R 

R 



n R V» , 

Pr V^V" - P R hf" 



(14) 
(15) 



Here, riR, pr, and Pr are the number density, the en- 
ergy density, and the pressure of the radiation fluid re- 
spectively. V is the velocity four- vector with V 2 = 1. 
h^ v = g^ v — V^V V is the projection operator on the 
plane perpendicular to V. The radiation fluid variables 
only depend on the temperature of radiation, Pr, since 
there are no relevant conserved quantum numbers besides 
R-parity, which is taken into account in the CDM fluid. 

The current density and the energy-momentum tensor 
of the imperfect CDM fluid can be written as 13 [l5 



J» = nU» + J (1) " , (16) 
T* v = p(T, n)U»U v - P(T, nW + P (1) . (17) 

n, U, p(T,n), and P(T,n) are the number density, four- 
velocity, energy density, and pressure of the CDM fluid, 
respectively. We omit the subscript x for the CDM com- 
ponent in Sees. IV, V, and VI, since the results of these 
sections hold true for more general forms of WIMP CDM. 
The projection h is orthogonal to U here. We do not in- 
troduce two different symbols in the following, because 
it is always clear from the context to which velocity h is 
referring. For the CDM fluid P and n are independent 
variables, since R-parity, i.e., the number of neutralinos, 



is conserved. P is not necessarily identical to Pr, al- 
though this is the case when both fluids are in thermal 
equilibrium. In the adiabatic limit all space-time gradi- 
ents are negligible, i.e., jW = 0, PW = 0, and the CDM 
fluid has the same temperature and the same velocity as 
the radiation fluid, U — V. 

For an imperfect fluid described by ( |l6| ) and ( |l7j ) num- 
ber density, energy density and velocity are not defined 
uniquely. To fix this ambiguity we define the number and 
energy density by 



n = UpJ* 



p{T,n) = U^U v T^ 



(18) 
(19) 



such that the hydrodynamic and thermodynamic defini- 
tions of n and p coincide. The velocity is fixed to be the 
comoving velocity of the CDM particles 



U> 



{J\J X 



J" 



(20) 



This choice of the velocity corresponds to the one of 
Eckar t [|13[ , and was applied to relativistic fluids by Wein- 
berg |15fl. An alternative would be the choice of Landau 
]l4| |, where the velocity is fixed to coincide with the ve- 
locity of the energy-momentum flow. In other words (|l8| ) 
requires J*- 1 -* to be perpendicular to U . In the same sense 
(|19| ) requires to project on the plane perpendicular 
to U. Condition ( p0[ ) means that has to vanish. 
With this choice 

= nU" , (21) 
T^ v = p(T, n)U»U v - P(T, n)h^ + TW " v . (22) 

The construction of P*- 1 -* from a first-order formalism is 
given in Jlj|. The starting point is that has to be ex- 
pressed in the equilibrium variables and their gradients. 
The basic observation is that the variation of the entropy 
per particle, a, along the adiabatic flow is generated by 
the change of P*- 1 ) in the same space-time direction, i.e., 



nT a 



TJ 71(1) P-V 



(23) 



where the dot denotes the hydrodynamic derivative, () = 
U^Q^. Equivalently we may write for the entropy cur- 
rent four- vector S» = naU^ + 1/P P A P (1) x,t (as can 
be easily seen in the comoving frame, S° = na = s is 
indeed the entropy density and TS l = P*- 1 -* 0i is the non- 
adiabatic contribution to the energy-momentum flow, 
which is the heat flow) 



T 2 s x^ 



(24) 



As a consequence only space-time derivatives of P and 
U can occur in TW in order to keep the rate of entropy 
production positive for all fluid configurations. 

The perturbed energy-momentum tensor may be ex- 
pressed in terms of the heat-flow vector 



) lt = T t „-TU ll 



(25) 
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and the traceless shear tensor 



(26) 



F {1) {uj, n, x) = A{u, x) + B^{u, i)n" + 



C /Jll/ {u,x)[n^n v + -h^ )+.... (32) 



With these abbreviations we write [|15j 

T (l) l»» = £ h H» IJX + v h VP h ™ Wpa 

X (h x »U u + h Xl/ U») Q x , 



(27) 



where £, r\ and x are the coefficients for bulk viscosity, 
shear viscosity and heat conduction. 

These parameters need to be calculated in the frame- 
work of a non-equilibrium theory. 



V. KINETIC DESCRIPTION OF CDM 

In this and the following section we generalize a 
method by Straumann to calculate the coefficients 
of transport for a species of massive particles that decou- 
ple kinetically. In Ref. Jl(| the problem of the decoupling 
of radiation quanta was treated. 

Let F(p, x) be the distribution function of neutralinos. 



It is clear from the kinetic equation (30) that we need to 
know how the functional L operates on In order to 

solve this problem we note that F is defined to be invari- 
ant under Lorentz-transformations. Let G x be the group 
of all Lorentz-transformations leaving V(x) invariant at 
every space-time point x, i.e. G x is the little group with 
respect to V. G x is isomorphic to the Lie-group 50(3). 
Since F is invariant under G x at every space-time point 
x, F^ is invariant and ( |32] ) is an expansion into irre- 
ducible polynomials with respect to G x . From equation 
( pOj ) it follows that the linear functional L is a scalar with 
respect to G x . Therefore it operates on the irreducible 
subspace spanned by the polynomials in (32) as a multi- 
ple of the identity. Thus we can write 



L[F W ] = -oj [k A + KiB^n"- 



(33) 



F(p,x) is normalized in such a way that F(p,x) d 3 pd 3 x where the Kj (j e {0, 1,2}) are functions of oj and 



gives the number of quanta in the volume d 3 x centered 
at the space-time point x and three-momentum within 
d 3 p. We assume that the neutralinos are close to thermal 
equilibrium and make the ansatz 



F = F (0) + F (1 



with 



where 



(2tt) 3 



exp 



F W < F (o) 5 ( 28 ) 



(29) 



only. Note that in the case of the model for the collision 
integral discussed above Kj = t -1 for all j E {1, 2, 3}. 

Next we derive expressions for A, and in terms 
of V, Tr, a and Kj using the kinetic equation (|30|). In 
order to do this we have to define a measure df2y on the 
two dimensional surface 5 = {p : p 2 = M 2 ,p° > and 
oj =const.}. d£ly is normalized such that 



1 

An 



dfl v = 1 



(34) 



± 1 



Here, Tr is the local temperature of the radiation fluid 
and a is the local pseudo chemical potential of the neu- 
tralino. To first order in the collision time, is a 

solution of the kinetic equation 



(p-0) F (0) =L[FW] 



(30) 



L is supposed to be a linear functional in f' 1 ', In the 
linear regime one often uses L^ 1 )] = -wr _1 F^ as a 
realistic model for the collision integral. 

In consideration of Eckarts approach to the hydrody- 
namics of imperfect fluids (ll| we introduce a four- vector 
perpendicular to V, 



The irreducible polynomials in ( p0|) are orthogonal with 
respect to dQv and are normalized as follows 

— f driven" = -- W" v , (35) 
4tt J s 3 

— h^ v /i Q/3 > - - h**" h af3 . (36) 
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Now it is possible to project out every tensor in the ex- 
pansion of L{F(^}], equation (32). Taking moments of 



( |30| ) and using (|32|) we obtain 



(31) 



A 



with oj = p ■ V such that F^> can be considered as a 
function of oj, n and x or equivalent as a function of the 
projection of p in the direction of V and perpendicular 
to it. 

Following ]22] ] we may now expand Fp\uj,n,x) into 
polynomials in n 



1 

_ 1 oj& 



Tr 3 

2 

1 r 



V 



— a 

OJ 



1 JR 
1 CJ$' 1 



-Qx 



To 



-a A 



(37) 

(38) 
(39) 
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Here, $'(o;/Tr — a) denotes the external derivative of 
F(°\ui,x). In calculating the integrals we replaced (p ■ 
d) T R with uj{V ■ d) T R and p (p ■ 9) V 13 with \p\ 2 n {n ■ 
d) V?. 

The coefficients ([}?]) and d3§| ) depend on the varia- 
tions of Tr and a along the adiabatic flow and the direc- 
tional derivative of a in the plane perpendicular to the 
adiabatic flow. In order to make sure that the rate of 
entropy production along the adiabatic flow is positive 
for all kinematical configurations these derivatives need 
to be proportional to space-time gradients of Tr and V. 
Using the adiabatic relations derived in Appendix B we 
find 



A 



B„. = 



1 ^ 

^0 Tr 
dP\ 
d Pj n 
1 CJ<A>' 



«1 Tr 



1 - 



K 1 


2\ 2 




m 2 \ 


1 
2 


Ld 


w 


~^ J 









dP 
dn 

h"Q> 



Cn,v — 



1 1 



7' 



(40) 



(41) 
(42) 



K 2 Tr 2 

with the enthalpy w = p + P 



A. Current density 

In kinetic theory the current density of neutralinos is 
given by 



(43) 

Considering our ansatz (|28|), we may write 

jk n _ jfe(o) ^ _|_ j fe (!) f* j (44) 

with the definitions 

Jfe (0) „ = /" ^ p P F (0) (Wj x) = n ky„ (45) 



poo 

4tt / dw(w 2 - Af 2 ) 1/2 w$ 
iiK 



(46) 



and 



jfc(i) /i = / *_RpP F <.i) ( W) nj j,.) = An y/« + jM g _ (47) 



Po 

An is generated by the coefficient A whereas Jdifi is gen- 
erated by B : 



,oo 

An = 4tt / <\lo(u 2 - M 2 ) 1/2 uj A , 
Jm 

~ r A^lo 2 - M 2 ) . 
3 Jm 



J diff 



(48) 
(49) 



Let us rewrite the above expressions with help of the 
following notation: 



f™(T a , a; x) = -An / du;(u, 2 - M 2 f 2 ^—- 

JM -LR K a 



(50) 



Note that the mass dimension of fa'^ is simply i + j. In 
terms of these functions we obtain 

An = 



ft 



(1,2) 



9P\ _ (3,0)1 (M) fdP 
dp) n h 3 +/ ° [dn 



/A 1 - _J_ 

diff " 3T R 



/( 3,0) _ (3,-1)" ^ 

n- 



V\, (51) 



(52) 



B. Energy-momentum tensor 

The energy-momentum tensor of neutralino CDM is 
given by 



d 3 p 
Po 



p^p v F(u>, n, x) . 



(53) 



Again, in consideration of our ansatz we may write 

rj^k iiv rpk(0) fif _j_ rpk(l) [iv (54) 

with the definitions 

d 3 p 



J Po 

p k yflyv _ pk hi lV ) 



P = 



P = 





POO 


An 


L 


An 




T 


JM 



rOO 

/ dc^ 2 - M 2 ) 3 / 2 $ , 
Jm 



and 



jrfc(l) \iV — rp(X) IU> _|_ rp(l) flV ^_ rp(l) flu _j_ 



(55) 
(56) 

(57) 



(58) 



The labels A,B,C indicate which tensor in the expansion 
(32j) gives rise to the extra contribution. From ([52]) and 
( p - @ we find 



T 



(1) u. V 



T, 



(1) ^ 



Ap U^ - AT /i^ 11 ' , (59) 
V { "h^ x Q\ , (60) 



1 

37k 



,(3,1) ,(3,0) w 

h h n 
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c 
with 



(61) 
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Ap 



(62) 



P (l,3) fdP\ _ (3,1)1 (1 , 2) / dP 

dp) n /o 3 +/o [dn 



V 



,x > 



AP 



/o 1 dp 



OP 



/•(5.-1) 1 , f (3,0) ( dP 

/o 3 + /o [dn 



V 



(63) 



A • 



VI. COEFFICIENTS OF TRANSPORT 

In the following we calculate the coefficients of bulk 
and shear viscosity and the coefficient of heat conduction 
for neutralino CDM starting from the kinetic description. 

In Sec. Ill we have introduced Eckart's approach to 
describe imperfect fluids fl3|| . The number density and 
the energy density of the CDM fluid coincide with the 
corresponding quantities in the adiabatic limit, see Eqs. 
@ and @, and the velocity of the CDM fluid is fixed 
via the particle current, see Eq. (p(i|). These definitions 
together with the required space-time symmetries and 
the second law of thermodynamics determine the most 
general structure of and see Jig ], as given in 

Eq. (§). 

To compare the kinetic description from the previous 
section with th e ap proach of Eckart it is necessary that 
the conditions (|l8|)-(po|) arc fulfilled. Instead we find for 
the kinetic description 



Vf,J kf " =n k + An, 

k fj,v 

h^Jx 



VpV v T K = p K + Ap 



J diS 



(64) 
(65) 
(66) 



Due to the non-equilibrium dynamics, the kinetic number 
density and energy density do not coincide with Eckart's 
definitions, (|64|) and ([35]) are in conflict with ( |l8| ) and 
(|l9|). Equation (|6^) shows the existence of a diffusion 
current in the plane perpendicular to V. As a conse- 
quence the current density four-vector does not point to 
the space-time direction that is required by the approach 
of Eckart (|oj). In the following we consider the tem- 
perature and the number density to be the independent 
thermodynamical variables. 

Let us first establish the link between the current in 
the kinetic and the hydrodynamic descriptions. As a first 
step we make a transformation of the velocity, such that 
the diffusion current vanishes, 



V ^ U - (n k + Any 1 J dies 



which allows us to write 



= (n k + An)U» = nU^ 



(67) 



(68) 



from comparison with ( |l8| ) . 

Let us now turn to the energy-momentum tensor. The 
transformation ( |6"7j ) with n — n k + An generates an extra 
contribution to the heat conduction since 



p k V^V v 



p k U t "U" 



pr_ 

n 



J diS u 



o(j- 



diff) 



(69) 



It remains to find the relation between the kinetic, Q6E 
and the hydrodynamic, (|l9|), definition of the energy den- 
sity. The point is that the definitions of temperature in 
the approach of Eckart and in the kinetic theory under 
consideration are different In kinetic theory there 

is a unique way to define temperature as the tempera- 
ture of the leptons and photons which stay in thermal 
equilibrium during and after the kinetic decoupling of 
the neutralinos. In the approach of Eckart the tempera- 
ture was chosen such that the energy density agrees with 
the one in the adiabatic limit. Thus it is clear that the 
difference in the definitions should be generated by Ap. 
Since we are only interested in effects linear in the col- 
lision time we may expand Ap in a first order Taylor 
expansion. Solving this expansion for the difference in 
the temperatures 



Tti = T 



dp 

df 



dp 

dn 



An - Ap 



(70) 



Let us now rewrite the energy-momentum tensor as cal- 
culated in the kinetic theory in terms of T, n and U : 

T k ^(T, n) = p(T, n)U p -U v - P(T, n)h» v 



A 



/; - 



+T, 



dP 

dn J \dp 

w 

B ' ' 

(1) nv 



Ap- AP 



n 



c 



(T,n) 



(71) 



This expression can be compared to (^7|) and the trans- 
port coefficients can be extracted. We express them in 
terms of the functions /« 



(1,3) 



dP 
dp 



p(3,l) 



(1,2) 



dP 
dp 



fr ( 


dP* 


dn 


+ 9 /o 


-i) 




2 




~ 3 



(72) 



dn ) 
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T X =- 3 



f(5 -1) 
J2 



ft 



(3,1) 



n 



+ 



(3,-1) 



(73) 



(74) 
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Instead of equating the collision integral to our expansion 
([33|) and solving for the unknown functions K a (u, x) [a S 
{0, 1, 2}) we give a qualitative correct estimate. This can 
be achieved by using the following model for the collision 
integral 



L[F«1 



-U) T 



(75) 



which corresponds to K a (u>, x) = r _1 for all (a € 
{0, 1,2}). This model reflects the linear dependence on 
and generates the variation of in the direction of 
the adiabatic motion through the rate of elastic scatter- 
ings. Furthermore since neutralinos are non-relativistic 
at kinetic decoupling we use the following approximations 



f(i,j) ~ f(i,j) ~ i\\(MTy-^M j nT 



for odd i and 



-(47,(^) 3/2 -p(.-f 



(76) 



(77) 



In this case Ap and AP depend linearly on An: Ap 
MAn ,AP 5/3TAn. At any time step At the varia- 
tion of the number density due to non-equilibrium pro- 
cesses is given by the number of collisions during At, i. e., 
An = fir. 

We find in 0(t) and up to order T/M 



(78) 



It is interesting to note that \ ~ at this order since 
the contribution of the transformation ([37|) to the energy 
momentum tensor cancels T^ K 

On the first sight it might be surprising that heat con- 
duction vanishes and bulk viscosity is nonvanishing. The 
mentioned cancellation between Tjp and J^m indicates 
that the only possible mechanism to transport heat in 
the neutralino fluid is convection. Since the neutralinos 
are very slow and very sparse heat can neither be radi- 
ated nor conducted. We decided to use the frame that is 
comoving with the neutralinos (Eckart's approach [fl3|| ), 
thus there is no heat conduction here. In a single fluid 
bulk viscosity goes along with internal degrees of free- 
doms or with particle production or decay. In our sit- 
uation the number of neutralinos is conserved and they 
do not have any internal degrees of freedom which can 
dissipate energy. Nevertheless, the bulk viscosity is non- 
zero, the reason is that we are dealing with two fluids and 
the bulk viscosity describes the energy dissipation to the 
radiation fluid. There should be a corresponding term 
for the radiation fluid, however we can neglect this term 
since 3> nT at kinetic decoupling. The authors of 
Ref. jllj have incorrectly assumed that \ 7^ and C = 
in their work. Let us note that our result (x = and 
C 7^ 0) holds in general for any kind of WIMPy CDM. 



VII. COLLISIONAL DAMPING OF ACOUSTIC 
PERTURBATIONS 

The viscosity coefficients and the coefficient for heat 
conduction enter in the decay rate of acoustic perturba- 
tions, which we will study now. Following Weinberg Jll| 
let us start with a static homogeneous fluid with 



U = (1,0) , p,P,n,T = const 



(79) 



This fluid should leave the adiabatic limit but stay close 
to thermal equilibrium. As a consequence small pertur- 
bations will occur with the space-time dependence 

S(p,P,n,T,k-U) = ^p (1 \P^\n^\T^\k-U^ 



exp 



(itut) exp (~ik • 



(80) 



Note that the perturbation of the zeroth component of 
U should vanish in order to guarantee the normalization 
condition. 



Inserting (80) in the conservation laws for the number 
density, energy density and momentum we get a system 
of three linear algebraic equations 



A((,r)) (6T,Sn,k-5tfj =0 , 



(81) 



where we used x = and the matrix ^4(£, rf} is given in 
|l5f . The dispersion relation is provided by the require- 
ment 

deU(C, rj) = , (82) 
which yields to first order in the collision time (k = \k\) 



Re w = kv s , Im u> 



-L[(,r)}k 2 = -r[(,r,]. (83) 



The square of the isentropic sound speed, v s , is given by 
T ( dP\ f dP\ n fdP\ 



p + P \dT J n \dp J n p + P\dnJ T 
and the characteristic length L for absorption reads 



(84) 



L[(,V] 



2p ~ 2 M s 



(85) 



Note that the length scale of collisional damping is pro- 
portional to the relaxation time. The authors of |ll]] 
assume instead that the characteristic scale for acoustic 
absorption is given by the collision time. This is correct 
only if acoustic perturbations are smeared out after a 
single contact with the heat bath. We proved already in 
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Sec. Ill that in the case under consideration a huge num- 
ber of contacts with the heat bath is needed to establish 
equilibrium. 

Since the parameters of the fluid are slowly varying 
during the cooling to Tkd, the amplitude of an acoustic 
perturbation behaves like a WKB solution. The damping 
of density perturbations is given by 



5(k) = exp — 



t(T kd ) 



r(t)dt 



exp 



3 Ti 



kd 



Cphys 



X, \ /I =Jkc 



H 



(86) 



We integrate over the time interval [0,ikd] during which 
the fluid is close to thermal equilibrium and CDM den- 
sity perturbations evolve like (damped) sound waves. In 
principle ikd is a function of k, however we find that for 
modes of interest we can take ik to be independent of k. 
This follows as Re(w)r = v s k p h ys T < 1 is easily fulfilled 
for the subhorizon scales (fc p hys/-ff )(21d) < V u s ~ 10 4 
for typical MSSM masses, including all modes of interest. 

From (^6|) we can read off a typical wavelength for 
collisional damping 



^d(Tkd) 



2tt 



^kd -Rff(Tkd) j 



(87) 



where «kd = \/ST\d/M^ and Rh = 1/H denotes the 
Hubble radius. We find l d ps 0.06(0.05)i? ff (T kd ) for 
= 100(150) GeV and M L = 200(250) GeV. 
Instead of characterizing acoustic perturbations by 
their wavelength or wavenumber fc p h ys ~ a ~ it 
is more convenient to work with a constant in time — 
the rest mass M of neutralinos within a sphere of radius 
27r/fc p i 1 y S : 



M = ~tt (27r/fc phys ) 3 ra^M* . 
Using the definition of M we can write (pq) as 



S(M) = exp - (M d /M) 3 



where the mass scale of damping, M d , is given by 



(89) 



2 4 7T 4 



(TkA 
5 V 10 7 \ M x 



rix (?kd ) R H (Tkd ) 



2.6 x 10 -s (! GeV )' M . 

(M^Tkd) 3 / 2 X ° 



(90) 



Exploring the parameter space of the MSSM we find typ- 
ically M D « 10" 9 M Q , cf. Fig. |. 
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FIG. 



100 125 

M x [GeV] 

4. The damping scale Md in solar masses as a func- 
tion of the neutralino mass My- for different values of the 
slepton Mass M~ L = 150, 200, 300, 400 GeV. The slepton mass 
increases from the bottom to the top. 



VIII. FREE STREAMING 

For temperatures T < T\ s the neutralinos are colli- 
sionlcss so that the viscosity coefficients vanish. Each 
neutralino moves along a geodesic in space-time. This 
geodesic motion of neutralinos provides a second damp- 
ing mechanism: free streaming [^3|. If the proper dis- 
tance 1{ S (t) which a neutralino can travel along a geodesic 
in time t is larger then then the proper wavelength 
A p hys = 27r/fc p h ys of a perturbation at t any structure 
will be wiped out since the neutralinos will propagate 
from an overdense region to an undcrdcnsc region. The 
proper distance of free streaming for a neutralino at time 
t e [iis,*eq] is given by 



k s {a) 



«ls 



vis In 
a \ a\ s 



R H (a) 



(91) 



(88) where a\ s denotes the expansion factor and v\ s 



y/ZT\ s /M^ the average neutralino velocity at last scat- 
tering. Exploring the parameter space of the MSSM we 
find at equality typically U s {a cq ) ~ 10 -8 i?#(a cq ), which 
" 4 pc today. This length scale con- 
10- 7 Af Q , cf. Fig. | 



corresponds to 5 x 10" 
tains a mass of Mf s 

Let us compare the scales of the two distinct damping 
mechanism at equality, the time when structures begin 
to grow (before equality CDM perturbations grow loga- 
rithmically). We find for a ^> a\ s 



ks 

Id 



10 m ( ^ 

\ais , 



2tt 



(92) 



which gives l{ s /ld ~ 6 (or M{ s /Md ~ 220) at equality. 
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FIG. 5. The free streaming mass Mf s at matter-radiation 
equality in solar masses as a function of the neu- 
tralino mass for different values of the slepton mass 



Mi 



150,200,300,400 GeV. The slepton mass increases 



from the bottom to the top. 



IX. CONCLUSIONS 

We have shown in this paper that kinetic decoupling 
of neutralino dark matter leads to collisional damping at 
the scale 10 _9 Mq. This scale could be larger for certain 
regions in the MSSM parameter space, e.g., when the 
neutralino mass and one of the slepton masses (probably 
the stau) are nearly degenerate. In that case our tree- 
level expressions become singular, and are not applicable. 
We have pointed out that it is important to distinguish 
between the collision time and the relaxation time of neu- 
tralino CDM. The corresponding temperatures differ by 
about an order of magnitude, which can lead to a differ- 
ence of several orders of magnitude in the corresponding 
mass scales of damping. 

The process of collisional damping has been described 
by imperfect fluids, and we calculated the transport coef- 
ficients from kinetic theory, by generalizing the method of 
Straumann Jl6| in order to include massive particles. We 
found that bulk viscosity can not be neglected, whereas 
heat conduction is negligible in the process of kinetic de- 
coupling of neutralinos. 

After kinetic decoupling free streaming starts to smear 
out remaining perturbations on scales below 10 -7 M Q by 
the time of equality. Both scales are quite close, which 
shows that both mechanisms have to be considered in 
the calculation of the resulting power spectra for cold 
dark matter. We will present the corresponding transfer 
functions and power spectra elsewhere |24[] . 

These damping mechanisms provide a sharp (exponen- 
tial) cut-off in the power spectrum of CDM objects. Such 
a cut-off sets the scale for the very first objects that form 
in hierarchical structure formation. Although this might 
be impossible to observe directly, it might have impli- 
cations on the substructure of galactic halos and on the 
structure of CDM in void regions, where some of the first 
CDM clouds might have a chance to survive. The cos- 



mological and astrophysical consequences of this cut-off 
will be investigated in a forthcoming publication [E4j . 
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APPENDIX A: CROSS SECTIONS 

This appendix contains the exact scattering ampli- 
tudes at tree level for elastic scattering and annihilation 
processes and simplified formulas for the related cross 
sections. 

Let us begin with the squared transition matrix ele- 
ment | T | 2 (summed over final and averaged over initial 
spins) for B + {L, L} — > B + {L, L} expressed as a func- 
tion of the usual Mandelstam variables. Our notation is 
as follows: 7/ (it) is the scattering amplitude which de- 
scribes the exchange of a left-handed slepton in the u 
channel and so on. For the squared terms we find 



|7i(u) | 2 =4( 5 6 L ) 4 

\T l {s) | 2 =A{gb L f 

\%(u) | 2 =4( 5 c L ) 4 

|T r ( S )| 2 =4( ffCi ) 4 



'u-M B 2 

JD 



u- M~ 2 
- - rr. 



s - M L A 



u - M f 



u — M- 2 
- M~ 2 - rr 



M L 



where we assumed for simplicity that the masses of the 
left- and right-handed sleptons are equal. 
For the different interference terms we find 



Mu)T r \s 



16<? 4 (b L c L ) 



16 5 4 (b L c L y 



2m 2 , Mf 2 - M-H 

^ B B 

u-M L 2 ) (s-M L 2 
2m 2 r M~ 2 ~ M B 2 t 

u B B _ 

s-M, 2 ' 
2 / m L Mg 



(u-Mf) (s-M L 2 ) 



(b LCL y 



u - M- 2 ) 
m L Mg \ 

2m j 2 

L B 



m,H 



(u-M L 2 )( 



.s - M L 2 
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Summing up the squared and interference terms yields 



T | 2 = 4 5 4 (6 L 4 



s — M^ 2 — m L 2 
. S ~ M f , 



u — M^ 2 — 7n L 2 
, U - M ? , 



(2m L M s y - 2Mj t 



8.g 4 (b L c L f 

' 2m L M h 
s-M £ * 



+ 



2m L M i} 



M L 



{2m L M s ) - 2m L 2 t 



M L 



)(• 



M L 



Since mi denotes the mass of a standard model fermion 
we use mi ~ 0. Furthermore the binos are nonrelativistic 
so that s « M^ + 2MgEL to a very good approximation. 
Using the mentioned simplifications we find 



256 (b[ 




and for the elastic cross section the simplified formula 

2 



24 



G ¥ M 2 



ii" 



M L 



M i, 



E L 2 



The squared transition matrix element (summed over 
final and averaged over initial spins) for B + B — * F + F 
may be found from | T | 2 for elastic scattering processes 
by making the following modifications: s — * u, t — * s 
and u — > t. Expanding in mp/Mg and in the Lorentz 
invariant relative velocity v up to second order yields 



nTda = dp 



w 



dn 



Since da must be a perfect differential 

*■(£).-» (* 

follows. For adiabatic motion 
= nTcj 



(Bl) 



(B2) 



(B3) 



ar/„ n \dT 



or, using the conservation law d(nU) = 
T 



' 1 % . ^ 



A • 



(B4) 



The Gibbs-Duhem relation gives the variation of the 
pseudochemical potential a = [x/T &s 



da 



dP w dT 



(B5) 



For adiabatic motion and using (B2) and (B4) this yields 



T a 



dP 

dn 



U' 



x ■ 



(B6) 



Using (B5) and the relativistic generalization of the Euler 
equation we find 



1 n 



(B7) 



Thus the variation of a in the plane perpendicular to the 
adiabatic flow is generated by the projection of the heat 
current on this plane. 



G F M? 



w 



M i, 



§0/ 



Mi 



{M § vf 



where v is given by (v/2) 2 = 1 — (2Mj 3 ) 2 /s. 



APPENDIX B: ADIABATIC RELATIONS 

In this appendix we show how T and space-time gra- 
dients of a are related to U x x and Q^. 

The second law of thermodynamics gives the variation 
in the entropy per particle a as 
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